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^ . Abstract. I examine the effect of trying to impose a Dirichlet boundary condition on a scalar field 

Tlj" ' by coupling it to a static background. The zero point - or Casimir - energy of the field diverges 

^-H , in the limit that the background forces the field to vanish. This divergence cannot be absorbed 

■ into a renormalization of the parameters of the theory. As a result, the Casimir energy of a surface 

on which a Dirichlet boundary condition is imposed, and other quantities like the surface tension, 
which are obtained by deforming the surface, depend on the physical cutoffs that characterize the 
coupling between the field and the matter on the surface. In contrast, the energy density away from 
the surface and forces between rigid surfaces are finite and independent of these complications. 
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INTRODUCTION 

This talk is offered to Joe Schechter on the occasion of his 65th birthday. Joe has always 
^ ! looked carefully at things that others took for granted, and often discovered new and 

^ I interesting physics as a result. Here I take a look at what boundary conditions do to 

^ the high energy behavior of interacting fields, and find some unexpected results. I hope 

Joe finds the story interesting and the argument convincing! The work described here 
is a by-product of a large collaboration including Eddie Farhi, Noah Graham, Vishesh 
Khemani, Markus Quandt, Marco Scandurra, Oliver Schroder, and Herbert Weigel. A 
summary of these results can be found in Ref. lU] and in a forthcoming publication. 10] 
Obviously boundary conditions are a very convenient idealization in field theory. Con- 
ducting boundary conditions give an excellent description of the behavior of electric and 
magnetic fields near good metals. Bag boundary conditions give a very useful charac- 
terization of the low modes of the quark field in hadrons. However physical materials 
cannot constrain arbitrarily high frequency components of a fluctuating quantum field, 
so the use of boundary conditions to model the effect of vacuum fluctuations, for exam- 
ple in the study of the Casimir effect, requires careful examination. Boundary conditions 
also appear in brane world scenarios for physics beyond the standard model and in lattice 
implementations of quantum field theories. 



point of view is that interactions between fields and matter are fundamental 
and that boundary conditions can only be substituted when they can be shown to yield 
the same physics. Surprisingly, there are important cases where they do not. The issue is 
the appearance of divergences or cutoff dependence in calculations of the total quantum 
zero point energy (the "Casimir energy") of a field subject to a boundary condition. 
If boundary conditions are placed on the field ab initio, methods have been developed 
that appear to yield a finite, cutoff independent energy.] 3] In contrast, we write down 
a quantum field theory describing the interaction of the fluctuating field with a static 
background, representing matter, and take a limit involving the shape of the background 
and the strength of the interaction that produces the desired boundary condition on a 
specified surface. Since the initial field theory is renormalizable, a finite, renormalized 
Casimir energy can be defined and computed for any non-singular background. The 
renormalized Casimir energy diverges in the boundary condition limit, indicating that the 
physical vacuum fluctuation energy depends in detail on the properties of the material 
that provides the physical ultraviolet cutoff. Note that these divergences have nothing 
to do with the standard infinities of quantum field theory, which were removed by 
renormalization. Instead they arise because the background necessary to enforce the 
boundary condition has too much strength at high frequencies. 

I don't want the point to get lost in complicated algebra, so I will discuss a case so 
simple that the necessary calculations are elementary. I will consider a scalar field in 
one dimension, obeying the Dirichlet boundary condition, (|) = 0, at one or two points. 
The limitation of this simple example is that the novel effects show up in the Casimir 
energy, but not in the Casimir force. In higher dimensions the effects are more dramatic 
and do effect measureable quantities. However the calculations are harder, so I will only 
report our results which are described in detail elsewhere. ||lil2] Of course this subject 
has been treated before, but not, as far as I know, in quite the way that I will describe. I 
will mention some of the other treatments later in my talk. 



DIRICHLET POINTS - A TOY MODEL IN ONE DIMENSION 

As a warm up, consider what it takes to enforce the boundary condition, (|)(0) =0 on 
an otherwise non-interacting and massless scalar quantum field in one dimension. This 
is the problem of the "Dirichlet point". The equation of motion for (|) coupled to some 
otherwise inert potential, a(x), is 

+ Xo(x){^ + n?^ = a)^(|) 

for a mode with energy CO. The coupling X is defined by normalizing o so that / o(jc) = 1 . 

To get all eigenmodes of (|) to vanish at x = it is necessary to take o to be "sharp" 
and X to be "strong": 

a(jc) 5(jc) SHARP 

X ^ oo STRONG 

In the sharp limit (|)' suffers a discontinuity at x = 0, A(|)'|o = X(|)(0). A careful study of 
the limit ?i — > oo shows that the boundary condition (|)(0) =0 emerges for all (O. A priori 



one would expect such a strong interaction to have a dramatic effect on the dynamics of 
(|), in particular on sum over zero point energies with the boundary condition compared 
to the energy without. 

First, suppose the boundary condition is imposed at the outset on all modes. This 
is the standard approach. The vacuum fluctuation energy is the sum over zero point 
energies of (|) minus the sum without the boundary condition, 

Ei = ^Y,(^^-^^o) (1) 

The situation is shown in Fig.[Tt. The solutions to the free field theory are (|) ~ s'mkx and 

^0(0) =0 0(-a) = O (j){a) = 



X 



FIGURE 1. Dirichlet points: a) one, and b) two at a separation of 2a 

^ ~ coskx for ^ > 0. With the boundary condition (|)(0) = the solutions are (|) ~ sinkx 
and (|) ~ sin \k\x, also with ^ > 0. Since the spectra are identical, (O and (Oq simply cancel 
and leave 

Ei=0 (2) 

The tilde is there to remind us that the boundary condition was imposed ab initio. The 
treatment in the standard texts is more careful than this, but in essence the same. Next, 
apply the same methods to the case of two Dirichlet points located at .x; = ±a (see 
Fig- eh)- The textbook result is 

Tl 

E2{a) = -— for m = (3) 

Eqs. © and @ are actually bizarre and unacceptible. First, remember that these are 
the total energies for each configuration relative to the vacuum - we did not drop any 
terms. To see the problem, consider first the limit a^ °°: 

lim£2(«) =0 = 2£i (4) 

a— >oo 

which is fine. It confirms that two widely separated Dirichlet points do not interact. Now 
consider the limit a ^ 0, 

lim E2 = —°° = E\ (5) 

which seems to say that the energy of a single Dirichlet point, the result of two that 
coalesce as a ^ 0, is infinite. Lastly, note that all of these results are suspicious because 
the massless scalar field theory suffers from infrared divergences in one dimension, 
leading us to expect log divergences where none have been found. 

Now let us consider the same problem from the perspective of an interaction of (j) with 
matter. We make the minimal model: we couple (|) to a non-dynamical scalar background 
field, o(x), with the (super) renormalizable interaction. 



L = -d^^d^^ - -m^(^^ - Xa(x)(^2 ^ c(e)a(jc) 



(6) 



Where c(£)o(x) is a counterterm and e is a cutoff, for example the fractional part of the 
dimension in dimensional regularization. We could elevate o to be a dynamical field by 
endowing it with an action of its own. However the core of the problem can be studied 
without this complication. L describes a renormalizable quantum field theory. In one 
dimension only one Feynman diagram - the tadpole - is divergent. This divergence 
is cancelled by the counterterm. As usual there is some scheme dependence involved in 
renormalization. We choose the "no tadpole scheme" where c(£) is chosen to completely 
cancel the tadpole graph so that (o) = 0. Any other scheme can be related to ours by 
a trivial shift in o. This theory has infrared divergences as m 0, so we keep m 7^ 
throughout. 

It is an straightforward to compute the renormalized energy of any configuration, 
£■[0], relative to the vacuum, a = 0, and to show that it is finite for piecewise continuous 
<5[x). The crucial question is what happens to £[0] when we try to take o to be sharp 
and \ to be strong. Do we reproduce the results obtained when the boundary condition 
is imposed ab initio or not? The sharp limit is benign in one dimension, so we can 
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FIGURE 2. Smooth but sharply peaked backgrounds approximating Dirichlet points: a) one, and b) two 
at a separation of 2a. In the "sharp" limit the backgrounds approach delta functions. 

take <3{x) to be a "spike" a(jc) = d{x) without difficulty. The renormalized energy of an 
isolated spike of strength X (see Fig.|2^) is found to be 

1 ..o dog (1 + 1)- 1 
E,M = ±- dt ) > (7) 

It is easy to see that this integral converges for any X and m^O, but diverges logarith- 
mically as m ^ 0. Likewise, the renormalized energy of two spikes at .x; = ±a is given 

by, 

1 flog(l + ^ + ^(l-e-4-))-X 
E2{a, X,m) = — dt ^ == ^ (8) 

Since these results are central to the rest of the discussion, I've appended a derivation 
(of eq. ©) at the end of the talk. 

Let's submit these results to the same examination as before: as a ^ 00 we find 

lim £'2(a,?i,m) = 2£'i(?i,m) (9) 

which is fine; and as a ^ 0, 

lim£'2(a,X,m) = £'i(2X,m) (10) 



which is much more satifactory than eq. Also, both Ei and E2 diverge logarithmi- 
cally as m 0, as expected. So the renormalized energy in a sharp background is finite 
and passes all the tests. 

However, the punchline is that both Ei and E2 diverge like —XlnX as °o. So the 
total renormalized Casimir energy diverges as the coupling constant, X, is taken strong 
enough to impose the boundary condition on all modes. X can be regarded as a cutoff, 
since modes with (O ^ X are not affected by the interaction. So in this example in a 
physical material the total, renormalized vacuum fluctuation energy is strongly cutoff 
dependent. Note that E\ and E2 are renormalized energies. There are no counterterms 
available to absorb this divergence. 

The two approaches disagree on the total energy, but they agree on the energy density 
and the force between two points. For any background (5{x), the energy density diverges 
only where o(x) 7^ 0. For a spike at a the energy density remains finite for x 7^ ±a in the 
limit X ^ 00, and the limiting form agrees with the density calculated using the boundary 
condition a priori.^ Likewise, the force between the two sharp sources agrees with the 
boundary condition calculation as ?i ^ 0°, 



In fact no measurement of the properties of the two points can detect the infinite energy 
stored locally on the boundary. This, however, is special to one dimension. Notice also 
that the vacuum fluctuation energy is negative (~ —XlnX as X ^ 0°). In a more realistic 
context this energy would be more than overwhelmed by the positive contributions to 
the energy coming from the curvature of a, which goes like |o'p, and the potential 
energy of a(x) which would involve higher powers of o{x) beginning with the mass 
term ^m^ / dx(5^{x). I left those terms out, not because they aren't present, but because 
I was trying to define an abstract problem - the "Dirichlet-Casimir" problem. Having 
failed, it is clear that the total energy depends not only on the vacuum fluctuation energy 
of (|), but also on the energy stored in the material represented by a. The two cannot 
be separated - no abstract "Casimir energy" can be defined for the Dirichlet boundary 
condition in one dimension. 

To summarize the results for a scalar field in one dimension: 

• The renormalized vacuum fluctuation energy is well defined and finite for any 
(piecewise continuous) background o(x). It differs from the energy calculated by 
assuming a boundary condition ab initio. The renormalized energy is material 
{ie. X), dependent, and diverges if X is taken to infinity to impose the boundary 
condition on all modes. 

• The change in the vacuum fluctuation energy with rigid displacement of the bound- 
aries is finite and cutoff independent and can be calculated by imposing the bound- 
ary condition ab initio. 




for 




(11) 



PHYSICAL EFFECTS IN OTHER DIMENSIONS 



The result of the previous section would be a merely a curiosity if it could not be mea- 
sured. In higher dimensions worse divergences occur and they too are always confined to 
the domain where the background fields are non-zero. Thus in any experiment in which 
the material (= the background fields) are unchanged, the cutoff dependence will cancel 
out. A case in point is the standard Casimir force between parallel, grounded, conducting 
plates. The force is measured by displacing the plates rigidly but not deforming them. 
The cutoff dependent terms remain unchanged as the plates are moved and the resulting 
force is finite. The result is the same whether you impose the boundary condition at the 
beginning or start with a smooth background and impose the boundary condition in a 
limiting process. 

The situation is completely different, however, if the material must be deformed to 
display the physical effect. The "Casimir pressure" on a sphere is the most interesting 
example . To measure this it is necessary to compare the total Casimir energy of a sphere 
of radius R with that of a sphere of radius R + dR. If there are cutoff dependent terms 
associated with the material, their contribution to the energy changes as the area changes, 
and they contribute to the pressure. Thus, calculations of Casimir pressures will differ 
between the case where boundary conditions are applied ab initio and where they are 
realized as a limit of the coupling to a background field. 

To examine these issues quantitatively we have studied the problem of the "Dirichlet 
sphere" - the boundary condition = imposed on a massive (or massless) scalar 
field in D dimensions. * D = 2 is the "Dirichlet circle" and D = 3 is the sphere. As in 
the one dimensional example, we replace the boundary condition by the coupling to a 
non-dynamical scalar background, a, with an interaction of the form 



where the . . . denote a finite series of counterterms necessary to renormalized the per- 
turbative divergences generated by the (|) — o interaction. For D = 2 and D = 3 only the 
counterterms shown are needed. We normalized the source so its integral over space is 
one. 



so the coupling constant, X, has dimension [mass]^^^^ 

We then compute the one-loop effective energy {ie. the Casimir energy), Ed{R, A, X), 
for a background o which is peaked at r = R and has a thickness A. We have studied 
GauBians and square barriers. The renormalized Casimir energy is finite for fixed A, and 
X. However it diverges in the sharp limit, A — ^ (except in one dimension, where we 
saw that the sharp limit was benign), so that not even the sharp limit exists for D > 2. If 



* Here I am describing work reported briefly in Ref. jT|. A more complete discussion is nearing comple- 
tion in Ref. 1 2] 

^ This corresponds to the physical situation that the surface gets thinner as R is increased. 



?io(r)4)2 + ci (£)o + c2(e)o2 + . . . , 




we keep A fixed, and take oo the Casimir energy diverges as well. The nature of the 
divergences depends on D. 

1 2 

• For the Dirichlet circle (D = 2), E2 diverges like ^ log A at fixed A as A ^ 0. Thus 
the tension (the two dimensional analog of pressure) diverges logarithmically as the 
thickness of the "circle" goes to zero. 

This is particularly nicely illustrated by examining the energy density, e(r), in a 
spherical shell between r and r + dr. The necessary formalism was developed in 
Ref. [4] In Fig. we plot e(r) for a GauBian background of fixed strength as 
a function of its width (denoted w in the figures). As w ^ the energy density 
diverges in a non-uniform way: at any fixed r ^ R, it approaches a limit. However 
the closer to R, the slower the convergence, and as a result, the total energy diverges 
as w ^ 0. The non-uniform behavior is clear in Fig.|3j). 

• For the Dirichlet sphere (D = 3), the leading divergence in £3 goes like 
Thus the pressure diverges as the thickness of the sphere goes to zero. 
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FIGURE 3. Casimir energy density in a GauBian background of width w approximating a Dirichlet 
circle, a) The energy density; b) a closeup of the behavior near r = 0.8/m. 

We expect that these divergences signal that the actual energy depends on the physical 
properties of the material. A represents its thickness and X plays the role of a frequency 
cutoff (eg. the plasma frequency, oOp) since modes with (O ^ ?i are unaffected by o. Since 
the other stresses to which materials are subject also depend on such properties of the 
material, we conclude that the Casimir pressure on a surface cannot be defined in a 
useful way that is independent of the particular dynamical properties of the material. 

Since the Casimir energy for any piecewise continuous background, o, has an expan- 
sion in Feynman diagrams, and since the divergence structure of Feynman diagrams is 
very well understood, it is interesting to look for the origin of the divergences in the dia- 
grams. The Casimir energy is simply the sum of all one-loop diagrams for (|) circulating 
in the background o. We find that the divergences never come from the loop integrations. 
Renormalization takes care of those. Instead the divergences come from the integration 
over the Fourier components of the external lines, o{p) . For D = 2 the divergence comes 
from the renormalized two-point function. Since the two-point function is easy to calcu- 
late, the divergence is easy to display. For D = 3 the two-point function also generates 



the leading divergence. However even the three-point function also diverges (logarith- 
mically) for D = 3. Since the loop integration for the three-point function is manifestly 
convergent, it provides a very clear demonstration that the divergences as A — > cannot 
be renormalized away by any counterterm available in the continuum theory. Instead 
they are physical manifestations of the cutoff dependence of the Casimir energy for the 
Dirichlet sphere. 

Our result disagrees with calculations which assume a boundary condition ab ini- 
tio.^ However those calculations handle divergences in an ad hoc manner and obtain 
suspicious results: for example the Casimir energy of a Dirichlet sphere is claimed to be 
finite in odd dimensions (D = 1,3,5,7...) and to diverge in even dimensions. 

The study of field theories in the presence of boundaries is an old subject. The 
divergences that occur near boundaries have been studied by many authors and with 
powerful tools. [7] In particular, it is well known that there are serious divergences near 
boundaries. Symanzik showed that all the divergences in the presense of boundaries (in 
certain theories) can be systematically cancelled by the introduction of new counterterms 
that live on the boundaries. However this begs our question of whether the divergences 
can be cancelled with only the counterterms available in an underlying interacting 
quantum field theory without boundaries. As far as I know, the implementation of 
boundary conditions as a limit of renormalizable, interacting quantum field theories is 
new. 

Although we have only completed our study for the Dirichlet boundary condition, we 
have examined conducting boundary conditions as well and expect our conclusions to 
apply to that case as well. 

Our interest in vacuum fluctuation energies goes far beyond the question I have 
discussed today. We are interested in whether quantum fluctuations can generate novel 
and complex behavior like solitons, in realistic quantum field theories. However, our 
investigations led us to a close examination of the nature of Casimir effects and to a 
result surprising enough to be worth telling Joe on the occasion of his birthday! 



This work has been done in collaboration with N. Graham, V. Khemani, M. Quandt, 
M. Scandurra, O. Schroder, and H. Weigel. The work I have reported is in large measure 
theirs. My work is supported in part by the U.S. Department of Energy (D.O.E.) under 
cooperative research agreement #DF-FC02-94ER40818. 



Here is an outline of the calculation of the vacuum fluctuation energy of an isolated 
"spike", eq. ©. We start with the universally accepted formula for the sum over zero 
point energies of a field, (|), fluctuating in a background o: 
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This is just the generalization to the continuum of j J^^oo — HqOq (with ^ = 1 of course). 
The (Oj are the energies of possible bound states. 6p(^) is the change in the density of 
states due to the background o. 5p is given by another well known result: 

where 6^(fc) is the phase shift for scattering in the partial wave in the background o 
(assumed to be symmetric). 

E[(5] defined by eq. (fT2b diverges in one dimension. To keep control of divergences we 
use dimensional regularization: we imagine that we are computing in D dimensions. It 
can be shown that eq. (fT2t is finite for < D < 1 . After isolating the possible divergences 
and renormalizing we will analytically continue to D = 1 . 

Next we use Levinson's theorem - which equates n times the number of bound states 
to the difference of the phase shift at /r = and k = oo (valid in D dimensions) to rewrite 
E[a] as 

^M = 2l(»^-^) + 2^y, dk{Vk^2-m)l^^ (13) 

j ^ 

At this point it is possible to show that the successive terms in the Born expansion of 5^ 
(the expansion of 5 in powers of o) can be put into one-to-one correspondence with the 
one loop Feynman diagrams.** In particular the first Bom approximation is identical to 
the tadpole diagram which diverges as D ^ 1 . 

To renormalize, we first subtract the first Bom approximation from eq. dT^ and add 
back the tadpole diagram to which it is equal. We then cancel the tadpole diagram 
against the counterterm c(£) (see eq. Q). The resulting renormalized and manifestly 
finite expression for the Casimir energy can now be evaluated dXD = 1, where the sum 
over phase shifts includes only the symmetric and antisymmetric channels: 

E[(5\ = ^ r dki^k^ + n?-m)^ ( h+ik)+h-ik)-h'^^\k)] (14) 
2n Jo dk \ / 

5*^^)(fc) is the Born approximation to the sum of the phase shifts. To save space I have 
dropped the sum over bound states since there are none in the case at hand. 

Now to the specific case of the 6-function background: The phase shift in the anti- 
symmetric channel vanishes. The symmetric channel phase shift is easily computed: 

5+(fc)=tan-i^ 

and the first Bom approximation is 

5"'W = 4 



** This is not true until after the Levinson's subtraction has been made. For a discussion and references, 
see H- 



The resulting integral, 



dk 



y/k^ + ni'^ — m 



k2{X^ + 4k^) 



is most easily evaluated by rotating the contour to the positive imaginary axis and 
integrating by parts. The result is the expression quoted in eq. Q. 
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